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We have completed a genome scan of a 12 -generation, 3,400-member pedigree with schizophrenia. Samples from 
210 individuals were collected from the pedigree. We performed an "affecteds-only" genome-scan analysis using 
43 members of the pedigree. The affected individuals included 29 patients with schizophrenia, 10 with schizoaffective 
disorders, and 4 with psychosis not otherwise specified. Two sets of white-European allele frequencies were 
used — one from a Swedish control population (46 unrelated individuals) and one from the pedigree (210 individuals). 
All analyses pointed to the same region: D6S264, located at 6q25.2, showed a maximum LOD score of 3.45 when 
allele frequencies in the Swedish control population were used, compared with a maximum LOD score of 2.59 
when the pedigree's allele frequencies were used. We analyzed additional markers in the 6q25 region and found a 
maximum LOD score of 6.6 with marker D6S253, as well as a 6-cM haplotype (markers D6S2S3-D6S264) that 
segregated, after 12 generations, with the majority of the affected individuals. Multipoint analysis was performed 
with the markers in the 6q25 region, and a maximum LOD score of 7.7 was obtained. To evaluate the significance 
of the genome scan, we simulated the complete analysis under the assumption of no linkage. The results showed 
that a LOD score >2.2 should be considered as suggestive of linkage, whereas a LOD score >3. 7 should be considered 
as significant. These results suggest that a common ancestral region was inherited by the affected individuals in 
this large pedigree. 



Introduction 

Mapping disease loci that predispose to mental disorders 
is difficult because of several obstacles, including non- 
Mendelian inheritance, phenotype definition, and gene 
and allele heterogeneity. One way to overcome these 
barriers is to study a single extended pedigree. The pre- 
sent study considers one of the world's largest reported 
pedigrees with individuals suffering from schizophrenia 
(MIM 181500). The pedigree has been traced back to 
the 1 7th century and comprises -3,400 individuals span- 
ning 12 generations. Here we report a complete genome 
scan for 43 affected individuals and 167 close relatives 
collected from the pedigree (in fig. 1, genotyped indi- 
viduals are indicated by gray boxes). Genotyping was 
done with 371 markers. Four different phenotypic cri- 
teria were used for the linkage analysis: (i) schizophrenia 
only (Scz); (ii) schizophrenia and schizoaffective de- 
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pressed type (SAD); (iii) schizophrenia, schizoaffective 
depressed type, and psychosis not otherwise specified 
(NOS) (Brdl); and (iv) all the previously mentioned di- 
agnoses and schizoaffective bipolar type (Brd2). Only 
the affected members of the pedigree were included in 
the LOD-score calculations. A caveat for our analysis is 
that, in a large pedigree such as that presented here, it 
is not possible to type all individuals from all generations 
to accurately determine the phase for all alleles. There- 
fore, it is important to estimate accurate allele frequen- 
cies for the linkage analysis (Boehnke 1991; Ott 1992; 
Freimer et al. 1993; Nechiporuk et al. 1993). To tackle 
this problem, we generated two sets of Swedish allele 
frequencies. The first set of allele frequencies was ob- 
tained from 46 unrelated individuals from the same ge- 
ographic area (i.e., northern Sweden) inhabited by the 
pedigree; this group was used because a control sample 
from northern Sweden would tend to allow reliable es- 
timation of the frequencies of common alleles present in 
the pedigree. The second set of allele frequencies used 
was derived from the pedigree, including all affected and 
unaffected individuals (210 individuals); this group is 
biased toward a false-negative result both because it in- 
cludes a large number (i.e., 43) of sick individuals and 
because the penetrance of the putative schizophrenia 
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Figure 1 Summary drawing of a large pedigree with schizophrenia, from northern Sweden. The total pedigree consists of 3,400 individuals, 
520 of whom are shown. The branches drawn were selected because they contain affected individuals who have given consent to be included 
in the genome scan and their closest relatives. The upper part of the pedigree, drawn in black, shows the first generations in the 17th century 
and the connections between the four branches drawn below. We have obtained blood samples from and have genotyped all the 210 individuals 
indicated by gray boxes. For simplicity, deceased individuals are not indicated. 
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gene is nor known. Our analysis revealed, at 6q25, a 
candidate schizophrenia-susceptibility locus, which was 
investigated by a fine-mapping strategy. 

Subjects, Material, and Methods 

Pedigree Construction and Clinical Assessment 

A pedigree consisting of 3,400 individuals spanning 
12 generations was constructed on the basis of church 
records, old psychiatric records, national registers, and 
information from a Ph.D. thesis (Sjogren 1935). All 43 
patients included in the genome scan received, before the 
start of the project, diagnoses within the "core spectrum 
of schizophrenia"— that is, Scz, schizoaffective psycho- 
sis, or psychosis not otherwise specified, according to 
DSM-III-R and/or the ICD-9 diagnostic systems. The 
original diagnosis and the number of previous hospi- 
talizations were obtained from the official Swedish in- 
patient psychiatric register. The patients were contacted 
again after 1994 and were asked to participate in the 
project. From 1994 to 2000, 34 of the 43 patients were 
interviewed again. Of these 34 patients, 26 were inter- 
viewed by means of the Swedish version of the Diag- 
nostic Interviews for Genetic Studies (DIGS) (Nurnber- 
ger et al. 1994); for the other 8 patients, a reliable and 
complete DIGS analysis was not possible, because of 
their psychiatric condition at the time of the interview. 
The remaining 9 of the 43 patients could not be inter- 
viewed again. Of these 9 patients, 6 died after the col- 
lection of blood samples, and 3 could not be assessed, 
because of either dementia (1 case), deafness (1 case), 
or stroke (1 case). For the 17 cases for which a complete 
DIGS analysis was not done after 1994, information 
from psychiatric unstructured interviews, psychiatric re- 
cords, clinicians, and interviews with first-degree rela- 
tives, by means of the Family Interview for Genetic Stud- 
ies (FIGS) (Nurnberger et al.- 1994), was considered 
sufficient for confirmation of the official clinical diag- 
nosis. The relatives of affected individuals were inter- 
viewed face-to-face and/or over the telephone. The phe- 
notypes of the relatives were considered as "unknown" 
for the statistical analysis. All clinical data obtained after 
1994, including those from DIGS, unstructured inter- 
views, psychiatric records, and FIGS, were added to the 
Swedish version of the OPCRIT 3.3 (Operational Cri- 
teria Checklist) (McGuffin et al 1991). After 1994, all 
43 affected individuals received DSM-IV and OPCRIT 
diagnosis, with the following distribution: Scz (29 cases); 
schizoaffective disorder, depressed type (6 cases), schi- 
zoaffective disorder, bipolar type (4 cases); and NOS (4 
cases). The OPCRIT-based diagnosis and the DSM- 
IV-based diagnosis were in full agreement. All diagnoses 
were established by two psychiatrists independently. The 
patients included 21 females with a mean age, at the 



time of the blood extraction, of 55 years (54.9 ± 15.5 
years) and a mean age at onset of about 27 years (26.6 
± 8.9years ). The 22 males had a mean age of -50 years 
(50.4 ± 15.2 years) and a mean age at onset of -24 
years (24.2 ± 7.3 years). 

Blood was collected from the 43 affected individuals 
and their closest relatives, resulting in a total of 210 
blood samples that were included in the genome scan 
(see fig.'l, where genotyped individuals are indicated by 
gray shading). The collection of all the blood samples 
was done by the Research Unit at the Department of 
Clinical Sciences, Division Psychiatry, at Umea Univer- 
sity in Sweden. 



Genotyping 

We used 371 polymorphic microsatellite markers. In- 
formation about many of the primers used for genotyp- 
ing and about sequences, heterozygosity, and allele sizes 
were obtained from The Genetic Location Database 
(Collins et al. 1996) and The Genome Database. The 
information for some additional, improved markers, as 
well as linkage-map information on the markers used 
was extracted from the databases at Axys Pharmaceu- 
ticals. The average marker spacing in the genome scan 
was -10 cM. There were eight gaps >30 cM; these gaps 
do not correspond to regions previously suggested to 
contain schizophrenia loci. The alleles for the markers 
were assigned by means of an approach for high- 
throughput genotyping (Hall et al. 1996). Multiple ali- 
quots of the DNA extracted from the 210 individuals, 
with a concentration of 4 ng/ml, were dried in multiple 
96-well replica plates. PCR mixing and the dispensing 
of mixes into the 96-well plates with the patients' DNA 
were performed by a Packard Robotic System Multi- 
probe 204 (Packard Instruments). Each PCR reaction 
was 20 pi, containing 2 pi of 10 mM dNTP, 2 ^tl of 10 
x PCR buffer, and 1 M l 50 of mM MgCI 2 . The DNA 
was subject to PCR in a MJ PTC-100 ® (MJ Research). 
Markers were amplified separately, and then 8-16 mark- 
ers were pooled according to panel-set category. Pooling 
was performed by a CRS/MultiPette Work Cell (CRS 
Robotics). The pooled samples were then loaded onto 
an ABI 373 or an ABI 377 (Applied Biosystems), de- 
pending on the resolution necessary for the markers. 
Thirty-three panel sets resulted in a total of 371 markers. 
The results were analyzed by GENOTYPER version 2.1 
(Applied Biosystems). An additional genome scan was 
performed, with the same 371 markers, on 46 healthy 
individuals originating from the same geographic area 
inhabited by the pedigree. This was done to obtain allele 
frequencies in the Swedish control population, to be used 
in the linkage analyses. 
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Figure 2 Genomewide simulation analysis. The upper graph 
shows the number of expected MLS values at or above a given MLS 
threshold, under the assumption of no linkage. This is equivalent to 
the average number of expected false-positive hits at or above each 
threshold, plotted on a logarithmic scale. Suggestive MLS values (1 
false-positive per genome scan) and significant MLS values (1 false- 
positive per 20 genome scans) are indicated by, respectively, the upper 
and lower dotted lines. The lower graph shows the proportion of 
genome scans exceeding each MLS value, plorted on a logarithmic 
scale. The dotted line represents the P = .05 level. 

Statistical Analysis 

Four different groups of phenotypic criteria were used 
in the linkage studies: (i) Scz, (ii) SAD, (iii) Brdl, and 
(iv) Brd2. Healthy members of the pedigree, as well as 
individuals with other psychiatric diagnoses, were de- 
noted "phenotype unknown." 

We performed a two-point linkage analysis, using the 
MLINK program from the LINKAGE package (Cot- 
tingham et al. 1993). We did not use allele frequencies 
from The Genome Database, for two reasons: first, The 
Genome Database frequencies are, in most cases, derived 
from a relatively small number of chromosomes, which 
might result in inaccurate frequencies that can give spu- 
rious results; second, the analyses of the Swedish control 
population and of the pedigree were done with the same 



set of tools, ensuring a correct identification of alleles in 
both groups. For the calculations, we used the same 
model that had been presented in our previous repli- 
cation work on chromosome 6p23 (Lindholm et al. 
1999); in brief, a model assuming a dominant trait was 
used. The sporadic risk was set to 0.4%, the mutation 
rate was 10~ 6 , and a maximum penetrance of .4 was 
used (Levinson et al. 1996). To avoid the problem of 
low penetrance of schizophrenia, only the affected mem- 
bers were included in the calculations. Linkage analyses 
were also performed with an additional set of allele fre- 
quencies. The second set of allele frequencies was ob- 
tained from the distribution of alleles in the pedigree 
analyzed in the genome scan. For saturation mapping, 
a third set of allele frequencies, from 46 healthy indi- 
viduals from northern Sweden, was used, depending on 
the availability of samples at the time for the analysis. 

To evaluate the significance of the results, we per- 
formed simulations of the complete genome-scan data. 
This analysis was performed at the PDC (Parallel Data 
Center) in Stockholm. We simulated the complete ge- 
nome-scan data, for the two sets of frequencies and four 
diagnostic groups (eight models), 100 times. For each 
simulated replicate of the genome scan, we analyzed all 
markers by means of the eight models. We therefore 
performed 37,000 simulations and 296,000 LOD-score 
calculations under the assumption of no linkage. We 
have reduced the number of alleles for each marker to 
five. For each marker, we noted both the maximum LOD 
score (MLS) under each of the eight models and the MLS 
of all eight models. For all genome-scan replicates, we 
then calculated, for the eight-model MLS, the average 
(and SD) number of LOD scores that exceeded each 
threshold (fig. 2 A) and the proportion of genome-scan 
replicates that had at least one LOD score that exceeded 
each threshold (fig. IB), A suggestive LOD-score thresh- 
old was considered as being the eight-model MLS value 
expected once in each genome scan; a significant LOD- 
score threshold was the average eight-model MLS value 
that would be expected to occur, by chance, once in 20 
genome scans. A similar analysis has been published else- 
where (Sawcer et al. 1997). 

We performed a three-point linkage analysis, using the 
FASTLINK program from the LINKAGE package (Cot- 
tingham et al. 1993). Seven markers were included in 
the calculations. The genetic distances were taken from 
The Genetic Location Database (Collins et al. 1996). 
The EHPLUS package (Zhao et al. 2000) was used to 
estimate the difference between the frequency of the 
6q25 haplotype in 43 affected individuals and that in 
46 unrelated controls. 

Results 

LOD-score data on the four different diagnostic groups 
described above (see the Subjects, Material, and Meth- 
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Table 1 

Markers with MLS Values >1.5, for Allele Frequencies in the 
Swedish Control Population and for Allele Frequencies in the 
Pedigree, for Four Diagnostic Categories 



Results for Allele 
Frequencies from 



Results for 
Allele Frequencies 



UnKUMUMJMt 


Control Population 


from Pedigree 


AINU IVIAKNCK 


MLS (0*) 


Diagnosis 


MLS (fl 3 ) 


Diagnosis 


D2S337 


2.60 (.15) 


Scz 


.00 (.35) 


Scz 


D2S1S1 


1.72 (.15) 


Scz 


.95 (.15) 


Scz 


3: 


2.14 (.20) 


Brdl 


.34 (.30) 


Brd2 




1.72 (.15) 


Brdl 


.01 (.30) 


Scz 


D3S1304 


2.72 (.10) 


Brdl 


.57 (.20) 


Brd2 




3.30 (.15) 


Brd2 


.30 (.25) 


SAD 


D3S1S66 


1.76 (.10) 


Brd2 


.53 (.15) 


Brdl 


D3S1S80 

5: 


1.55 (.15) 


Scz 


.18 (.30) 


Scz 


DSS419 


2.02 (.20) 


Brdl 


.47 (.20) 


Brd2 


DSS644 


1.82 (.15) 


Scz 


.49 (.25) 


SAD 


D5S433 

6: 


1.59 (.20) 


Brd2 


.08 (.30) 


Scz 


D6S344 


3.60 (.05) 


Scz 


-.03 (.40) 


Scz 


D6S264 

7: 


3.45 (.00) 


Brdl 


2.59 (.00) 


Brdl 


D7S636 

8: 


2.00 (.30) 


Scz 


.00 (.20) 


Scz 


D8SSS0 

9: 


2.90 (.15) 


Scz 


.70 (.15) 


Scz 


D9S1S7 


1./2 (.05) 


iCZ 


1.47 (.05) 


Scz 


D9S17S 


2.00 (.15) 


SAD 


.00 (.20) 


SAD 


D9S299 


2.01 (.05) 


Brd2 


1.11 (.10) 


Brd2 


10: 










0 J OS J 89 
11: 
DI1S1324 


1.94 (.10) 


Brd2 


.38 (.20) 


Brd2 


174 (.20) 


Brd2 


-.03 (.40) 


Scz 


13: 










DJ35253 


1.66 (.15) 


Scz 


.65 (.15) 


Scz 


D J 35 J 52 


1.58 (.20) 


Brd2 


.04 (.35) 


Brdl 


D J 351 73 


1.53 (.10) 


SAD 


.53 (.15) 


SAD 


14: 










D14S25S 


1.60 (.15) 


Brdl 


.00 (.20) 


Brdl 


16: 










D16S40I 
20: 
D20S13 8 


2.52 (.15) 


SAD 


.23 (.30) 


Scz 


3.57 (.20) 


Brd2 


-.12 (.40) 


Scz 


D205175 


2.70 (.15) 


Brd2 


.10 (.25) 


Brdl 


D205100 


1.76 (.10) 


Scz 


1.37 (.10) 


Scz 


D20S173 
21: 
D215266 


2.39 (.15) 


SAD 


.05 (.35) 


Scz 


1.61 (.15) 


Scz 


.27 (.25) 


Scz 


22: 










D22S2S0 


1.61 (.15) 


Brd2 


.09 (.30) 


Brd2 



Note. — The MLS values for marker D6S264, which reached sug- 
gestive levels of significance in both the Swedish control population 
and the pedigree, are underlined. 

' Recombination fraction. 



ods section), analyzed by means of allele frequencies in 
the Swedish control population and in the pedigree, are 
plotted in figure 3; this figure shows that linkage analyses 
using allele frequencies in the Swedish control popula- 
tion gave LOD scores >3 with four markers: D3S1263, 
D6S344, D6S264, and D20S118. Of these, only 
D6S264 gave LOD scores >2 when allele frequencies in 
the pedigree were used (fig. 3B and table 1). To evaluate 
the significance of the results, we performed a complete 
simulation of the genome scan, using all models to cal- 
culate empirical LOD scores for the pedigree, under the 
assumption of no linkage. We have reduced the number 
of alleles for each marker to five. We performed 37,000 
simulations (370 markers simulated 100 times each) and 
296,000 LOD-score calculations (each of the simulated 
pedigrees was subjected to eight LOD-score calculations, 
for all combinations of allele frequencies and diagnostic 
groups). These simulations indicated that, for the ped- 
igree structure and models used in the present study, a 
LOD score >2.2 should be considered as suggestive 
whereas a LOD score >3.7 should be considered as sig- 
nificant (fig. 2). Table 1 shows that nine markers in- 
cluded in the genome scan produced LOD scores higher 
than the suggestive value of 2.2: markers D2S337, 
D3S1304, D3S1263, D6S344, D6S264, D8S550, 
D16S401, D20S118, and D20S173. Table 1 also in- 
cludes additional markers with LOD scores >1.5, for 
comparisons with other studies. All nine markers that 
produced suggestive linkage results were selected for fur- 
ther studies. Flanking markers on each side of these 
markers were typed. The results supported only the pos- 
itive linkage for D6S264. Two-point LOD scores for six 
additional markers located close to D6S264 are shown 
in figure 4. A two-point MLS of 6.6 was obtained with 
D6S253 (marked by an arrow in fig. 4). Figure 4 also 
shows haplotypes in the 6q25 region, constructed for all 
affected individuals whose parents were available. For 
19 affected individuals connected to the same ancestor 
. nine generations back, we were able to determine the 
phase and the haplotype inherited from the pedigree; 11 
of these 19 patients shared a 6-cM haplotype 
(D6S2S3-D6S264), and 7 of the 19 patients shared a 
part of this haplotype. 

We also conducted three-point analysis of the seven 
markers in the 6q25 region, using allele frequencies in 
the Swedish control population, and we found an MLS 
of 7.7 between markers D6S2S3 and D6S297 (fig. 5). 

Discussion 

We have conducted a genomewide screening of a very 
large schizophrenia pedigree from northern Sweden, and 
we have found, at 6q25, a candidate schizophrenia-sus- 
ceptibility locus. To evaluate the significance of the re- 
sults, the calculations should be corrected for the fact 



Lindholm et al.: Schizophrenia-Susceptibility Locus at 6q25 

^ LOD scores obtained with Swedish allele frequencies 



BEST AVAILABLE COPY °' 




.4; w 

•X ; ; ; 


Diagnosis 


;t- : : ". 


•: ♦ ;' : 

: & ; 


♦ MaxLODScz 




♦ Max LOD SAD 




▲ Max LOD Brd1 




xMaxL0DBrd2 



B 



LOD scores obtained with pedigree allele frequencies 



-2 





D6S264 










1 1 


'■ '. ': ! 






▲ 








! • 


f s : -: 






: ; ** x 


. i I 




♦ TS 


» \* 

: / : : 





2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 202U2 X 

Chromosome 



Diagnosis 
*MaxLODScz 
♦ Max LOD SAD 
a Max LOD Brdl 
x Max LOD Brd2 



Figure 3 Plot of genome-scan data on 371 markers, using four diagnostic criteria and two sets of allele frequencies. The MLS for each 
marker is plotted for allele frequencies in the Swedish control population [A) and in the pedigree (J3). Each panel includes one plot for each of 
the four diagnostic categories used (four plots per panel): Scz, SAD, Brdl, and Brd2. The plotted values for the LOD scores >3 are indicated 
by an arrow and the marker name. The scale for each chromosome reflects the number of markers typed, rather than the size of the chromosome. 



that we performed multiple tests. A very stringent cor- 
rection would be to decrease the value of all the LOD 
scores by log(8) = 0.9, since we have used four types 
of diagnosis and two sets of allele frequencies. This cor- 
rection, which would render all of our reported LOD 
scores as not significant, is not very appropriate, since 
the four phenotypes are correlated. As an alternative, 
we performed a complete simulation of the genome scan, 
using all eight models to calculate empirical LOD scores 
for the pedigree, under the assumption of no linkage. 
These simulations indicated that, with the pedigree struc- 
ture and models used in the present study, a LOD score 
>2.2 should be considered as suggestive whereas a LOD 
score >3.7 should be considered as significant. As men- 
tioned above, the highest linkage value in the genome 



scan was reached with marker D6S264 on chromosome 
6q25 (table 1). The results for this marker reached sug- 
gestive values of significance, according to our simula- 
tion studies. The fact that D6S264 was the only marker 
that gave consistent results with the two sets of allele 
frequencies further strengthens the significance of these 
results. 

Additionally, D3S1263 located on 3p25, D6S344 lo- 
cated on 6p24, and D20S118 located on 20pll.2 
showed the highest suggestive values (MLSs of 3.30, 
3.60 and 3.57, respectively). MarkerD653 44 gave an 
MLS of 3.6. This marker is located on 6p24, a region 
that others (Straub et al. 1995; Schizophrenia Linkage 
Collaborative Group for Chromosomes 3, 6 and 8 
1996; Turecki et al. 1997; Nurnberger and Foroud 
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Figure 4 Haplotype analysis of the 6q25 region. In the table, a haplorype (indicated by gray shading) close to marker D6S264 is shown, 
for 19 affected individuals whose phase is known. Only the haptotype that is transferred from the pedigree (and not from the in-married parent) 
is shown. For the individuals in whom a marker was uninformative, the allele included in the haplotype is indicated by stippling. Both alleles 
are indicated for genorypes for which phase is unknown. Haplorypes from siblings are not separated by a vertical black line. The tale also 
shows, for each marker, the LOD score, the recombination fraction (0), and the diagnosis that gave the MLS. The circles to the right of the 
table show allele frequencies in affected individuals and in the Swedish control population. Alleles that are represented in the haplotype are 
indicated by gray shading. 



1999) have suggested as containing a schizophrenia- 
susceptibility locus. In fact, we elsewhere have published 
a part of the 3,400-member pedigree (i.e., -2,000 mem- 
bers), in an attempt to replicate findings for chromo- 
some 6p (Lindholm et al. 1999). D6S344 is located -2 
cM from marker D6S277, contained in a haplotype that 
was found to segregate with the disease in a small 
branch of the pedigree (Lindholm et al. 1999); however, 
one possibility is that D3S1263 (LOD score 3.3), 
D6S344 (LOD score 3.6), and D20SU8 (LOD score 
3.6) are false-positive results, since the calculations were 
not consistent for both sets of allele frequencies. The 
results stress the importance of an accurate evaluation 
of allele frequencies used in the analysis of extended 
pedigrees. Another possibility is that there is hetero- 
geneity inside the pedigree, with one or more of these 
markers segregating with a schizophrenia locus in just 
one part of the pedigree. 

LOD scores were calculated for additional markers 
in the 6q25 region, and several of them reached sug- 
gestive or significant values, as shown in figure 4. One 
problem in this region is that the alleles present in the 
haplotype segregating in the pedigree are frequent in the 
Swedish control population, a finding that decreases 
their informativeness for LOD-score calculations (in fig. 
4, circles represent allele frequencies in the affected in- 
dividuals in the pedigree and in a Swedish control pop- 



ulation). However, marker D6S253 (indicated by an 
arrow in fig. 4), which has a more-even allele distri- 
bution, resulted in an MLS of 6.6. According to the 
three-point analysis (fig. 5), it is possible that a schizo- 
phrenia-susceptibility gene is located in the region be- 
tween markers D6S253 and D6S297, a region that 
reached an MLS of 7.7. The haplotypes shown in figure 
4 supports this finding, except in the case of patient 
135; it is possible that allele 2 for marker D6S297 has 
mutated to allele 3 in this individual; an alternative 
explanation is that this individual is a phenocopy. All 
other individuals have either the complete haplotype ( 1 1 
of 19 patients for whom we could determine haploty- 
pes) or part of it (7 of 19 patients). We were not able 
to collect, for haplotype analysis, all the parents of the 
43 affected individuals. Therefore, we used the EHPLUS 
package (Zhao et al. 2000) to estimate the frequency 
of the D6S253 (allele 2)-D6S297 (allele 2)-D6S264 
(allele 4) haplotype in the 43 affected individuals from 
the pedigree versus that in 46 unrelated controls. We 
found a significant increase in the frequency of this 6- 
cM haplotype in the pedigree compared with the con- 
trols: empirical P = .004 for a model-free analysis; in 
light of the frequency of the haplotype in the normal 
population, this P value indicates a significant proba- 
bility that this region has been inherited by the affected 
individuals IBD. This result suggests that the increase 
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Figure 5 Multipoint analysis of the markers in the 6q25 region, 
performed- by use of diagnostic category Brdl and allele frequencies 
in the Swedish control population. Three-point LOD scores were cal- 
culated and plorted against the distance (in cM) from the p-terminal. 
The MLS was obtained berween markers D6S2S3 and D6S297 and 
reached a value of 7.7. 

in haplotype frequency in the pedigree is not due to a 
common frequency of this region in the Swedish pop- 
ulation. In conclusion, the EHPLUS analysis further 
strengthens the possibility that this region segregates 
with the disease in the large pedigree. 

Previous schizophrenia genome scans have proposed 
schizophrenia-susceptibility loci on lq, 2q, 3p, 4q, 5p, 
5q, 6p, 7q, 8p, 9q, lOp, lOq, llq, 13q, 22q 5 and Xp 
(Lasseter et al. 1995; Pulver et al. 1995; Straub et al. 
1995, 1997, 1998; Gill et al. 1996; Blouin et al. 1998; 
Crowe and Vieland 1998; Levinson et al. 1998; Brzus- 
towicz et al. 1999, 2000; Hovatta et al. 1999; Nurn- 
berger and Foroud 1999; Ekelund et al. 2000). Two of 
these genome scans reported significant linkage results 
(Blouin et al. 1998; Brzustowicz et al. 2000). These two 
studies used dominant and recessive models for the seg- 
regation analysis. In the present study, we used only a 
dominant model, because we considered a dominant 
segregation to be more likely in the large pedigree. The 
analysis of only a dominant model decreases the risk of 
type 1 errors. The studies by Blouin et al. (1998) and 
Brzustowicz et al. (2000) included the analysis of mul- 
tiplex families, whereas the present study includes a sin- 
gle large pedigree. The analysis of a single pedigree de- 
creases the probability of heterogeneity and increases 
the chance of finding a single gene of major effect. Mul- 
tiple susceptibility loci are most probably involved in 
the etiology of schizophrenia, and the fact that these 
three studies report different loci is not surprising. 

Until now, a schizophrenia locus at 6q25 had not been 
described; however, a 6q position closer to the centro- 
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mere has been proposed to contain a schizophrenia- 
susceptibility locus (Cao et al. 1997; Kaufmann et al. 
1998; Martinez et al. 1999). In fact, excess allele sharing 
has been reported for markers in the 6q 13-26 region 
(Cao et al. 1997; Kaufmann et al. 1998); however, the 
6q23-q26 region was excluded in a follow-up multi- 
center study, and only 12 microsatellite markers span- 
ning 6ql3-q23 were investigated (Martinez et al. 1999). 
The 6q23-q26 region was most likely excluded because 
the P values for the markers in this region were lower 
than those for the markers on 6ql3-q23. The possibility 
remains that all the studies of the 6q region describe a 
single schizophrenia locus. 

Interestingly, an autism-susceptibility locus also has 
been mapped to the 6q25 region (Philippe et al. 1999), 
and deletions of the 6q25 segment can result in devel- 
opmental problems and anomalies of the brain (Suku- 
mar et al 1999). These results suggest that genes im- 
portant for normal brain function and development are 
located on 6q25. 

In summary, we have analyzed the world's largest 
reported pedigree with individuals affected with schizo- 
phrenia and have found a candidate schizophrenia-sus- 
ceptibility locus at 6q25. The fact that the 6q locus 
reached one of the highest LOD scores reported so far 
warrants further investigations of this region, in other 
families around the world. 

Acknowledgments 

This work was supported by funding from Axys Pharma- 
ceuticals and Parke-Davis (1997-2000); by Swedish Medical 
Research Council grants 21X-10412, KO0-Z1X-13443-01A, 
and K1999-27X-012695-02B; by theTorsten and Ragnar Sod- 
erbergs Foundation; and by the Beijer Foundation. We are 
grateful to Dr. Joe Terwillinger, for his help in the design of 
the genomewide simulation analysis, and to Per Oster and the 
rest of the staff at the parallel computer center in Stockholm, 
for their invaluable help with the calculations. We would like 
to acknowledge Betsy Nanthakumar, who organized the use 
of the genotyping platform. We also wish to thank the patients 
and their families who participated in this study. E.J is partially 
supported by a grant from the Goran Gustavsson Foundation. 



Electronic-Database Information 

The accession number and URLs for data in this article are 
as follows: 

Genetic Location Database, The, http://cedar.genetics.soton.ac 

.uk/public_html/ldb.html (for primers used for genotyping, 

sequences, heterozygosity, and allele sizes) 
Genome Database, The, http://www.gdb.org/ (for primers used 

for genotyping, sequences, heterozygosity, and allele sizes) 
Online Mendelian Inheritance in Man (OMLM), http://www 

.ncbi.nlm.nih.gov/Omim/ (for schizophrenia |MIM 181500]) 



104 



Am.]. Hum. Genet. 69:96-105, 2001 



PDC (Parallel Data Center [Center for Parallel Computers, 
KTH]), http://www.pdc.kth.se 

References 

Blouin JL, Dombroski BA, Nath SK, Lasseter VK, Wolyniec 
PS, Nestadt G, Thornquist M, et al (1998) Schizophrenia 
susceptibility loci on chromosomes 13q32 and 8p21. Nat 
Genet 20:70-73 

Boehnke M (1991) Allele frequency estimation from data on 
relatives. Am J Hum Genet 48:22-25 

Brzustowicz LM, Hodgkinson KA, Chow EW, Honer WG, 
Bassett AS (2000) Location of a major susceptibility locus 
for familial schizophrenia on chromosome Iq21-q22. Sci- 
ence 288:678-682 

Brzustowicz LM, Honer WG, Chow EW, Little D, Hogan J, 
Hodgkinson K, Bassett AS (1999) Linkage of familial schizo- 
phrenia to chromosome 13q32. Am J Hum Genet 65: 
1096-1103 

Cao Q, Martinez M, Zhang J, Sanders AR, Badner JA, Crav- 
chik A, Markey CJ, Beshah E, Guroff JJ, Maxwell ME, 
Kazuba DM, Whiten R, Goldin LR, Gershon ES, Gejman 
PV (1997) Suggestive evidence for a schizophrenia suscep- 
tibility locus on chromosome 6q and a confirmation in an 
independent series of pedigrees. Genomics 43:1-8 

Collins A, Frezal J, Teague J, Morton NE (1996) A metric map 
of humans: 23,500 loci in 850 bands. Proc Natl Acad Sci 
USA 93:14771-14775 

Cottingham RW Jr, Idury RM, Schaffer AA (1993) Faster se- 
quential genetic linkage computations. Am J Hum Genet 53: 
252-263 

Crowe RR, Vieland V (1998) Chromosome 5 workshop. Psy- 
chiatr Genet 8:73-78 

Ekelund J, Lichtermann D, Hovatta I, Ellonen P, Suvisaari J, 
Terwilliger JD, Juvonen H, Varilo T, Arajarvi R, Kokko- 
Sahin ML, Lonnqvist J, Peltonen L (2000) Genome-wide 
scan for schizophrenia in the Finnish population: evidence 
for a locus on chromosome 7q22. Hum Mol Genet 9: 
1049-1057 

Freimer NB, Sandkuijl LA, Blower SM (1993) Incorrect spec- 
ification of marker allele frequencies: effects on linkage anal- 
ysis. Am J Hum Genet 52:1102-1110 

Gill M, Vallada H, Collier D, Sham P, Holmans P, Murray R, 
McGuffin P, et al (1996) A combined analysis of D22S278 
marker alleles in affected sib-pairs: support for a suscepti- 
bility locus for schizophrenia at chromosome 22ql2. Schizo- 
phrenia Collaborative Linkage Group (Chromosome 22). 
Am J Med Genet 67:40-45 

Hall JM, LeDuc CA, Watson AR, Roter AH (1996) An ap- 
proach to high-throughput genotyping. Genome Res 6: 
781-790 

Hovatta I, Varilo T, Suvisaari J, Terwilliger JD, Ollikainen V, 
Arajarvi R, Juvonen H, Kokko-Sahin M-L, L Vaisanen, 
Mannila H, Lonnqvist J, Peltonen L (1999) A genomewide 
screen for schizophrenia genes in an . isolated Finnish sub- 
population, suggesting multiple susceptibility loci. Am J 
Hum Genet 65:1114-1124 

Kaufmann CA, Suarez B, Malaspina D, Pepple J, Svrakic D, 
Markel PD, Meyer J, Zambuto CT, Schmitt K, Matise TC, 
Harkavy Friedman JM, Hampe C, Lee H, Shore D, Wynne 



D, Faraone SV, Tsuang MT, Cloninger CR (1998) NIMH 
Genetics Initiative Millenium Schizophrenia Consortium: 
linkage analysis of African-American pedigrees. Am J Med 
Genet 81:282-289 

Lasseter VK, Pulver AE, Wolyniec PS, Nestadt G, Meyers D, 
Karayiorgou M, Housman D, et al (1995) Follow-up report 
of potential linkage for schizophrenia on chromosome 22q: 
part 3. Am J Med Genet 60:172-173 

Levinson DF, Mahtani MM, Nancarrow DJ, Brown DM, 
Kruglyak L, Kirby A, Hayward NK, Crowe RR, Andreasen 
NC, Black DW, Silverman JM, Endicott J, Sharpe L, Mohs 
RC, Siever LJ, Walters MK, Lennon DP, Jones HL, Nertney 
DA, Daly MJ, Gladis M, Mowry BJ (1998) Genome scan 
of schizophrenia. Am J Psychiatry 155: 74 1-750 

Levinson DF, Mowry BJ, Sharpe L, Endicott J (1996) Pene- 
trance of schizophrenia-related disorders in multiplex fam- 
ilies after correction for ascertainment. Genet Epidemiol 13: 
11-21 

Lindholm E, Ekholm B, Balciuniene J, Johansson G, Castens- 
son A, Koisti M, Nylander P, Pettersson U, Adolfsson R, 
Jazin E (1999) Linkage analysis of a large Swedish kindred 
provides further support for a susceptibility locus for schizo- 

. phrenia on chromosome 6p23. Am J Med Genet 88: 
369-377 

Martinez M, Goldin LR, Cao Q, Zhang J, Sanders AR, Nan- 
carrow DJ, Taylor JM, Levinson DF, Kirby A, Crowe RR, 
Andreasen NC, Black DW, Silverman JM, Lennon DP, Nert- 
ney DA, Brown DM, Mowry BJ, Gershon ES, Gejman PV 
(1999) Follow-up study on a susceptibility locus for schizo- 
phrenia on chromosome 6q. Am J Med Genet 88:337-343 

McGuffin P, Farmer A, Harvey I (1991) A polydiagnostic ap- 
plication of operational criteria in studies of psychotic ill- 
ness: development and reliability of the OPCRIT system. 
Arch Gen Psychiatry 48:764-770 

Nechiporuk A, Fain P, Kort E, Nee LE, Frommelt E, Polinsky 
RJ, Korenberg JR, Pulst SM (1993) Linkage of familial Alz- 
heimer disease to chromosome 14 in two large early-onset 
pedigrees: effects of marker allele frequencies on lod scores. 
Am J Med Genet 48:63-66 

Nurnberger JI Jr, Blehar MC, Kaufmann CA, York-Cooler C, 
Simpson SG, Harkavy-Friedman J, Severe JB, Malaspina D, 
Reich T (1994) Diagnostic interview for genetic studies: ra- 
tionale, unique features, and training. NIMH Genetics Ini- 
tiative. Arch Gen Psychiatry 51:849-859; discussion 
863-864 

Nurnberger JI Jr, Foroud T (1999) Chromosome 6 Workshop 

report. Am J Med Genet 88:233-238 
Ott J (1992) Strategies for characterizing highly polymorphic 

markers in human gene mapping. Am J Hum Genet 51: 

283-290 

Philippe A, Martinez M, Guilloud-Bataille M, Gillberg C, Ras- 
tam M, Sponheim E, Coleman M, Zappella M, Aschauer 
H, van Malldergerme L, Penet C, Feingold J, Brice A, Le- 
boyer M (1999) Genome-wide scan for autism susceptibility 
genes: Paris Autism Research International Sibpair Study. 
Hum Mol Genet 8:805-812 

Pulver AE, Lasseter VK, Kasch L, Wolyniec P, Nestadt G, 
Blouin JL, Kimberland M, Babb R, Vourlis S, Chen H, Lal- 
ioti M, Morris MA, Karayiorgou M, Ott J, Meyers D, An- 
tonarakis SE, Housman D, Kazazian HH (1995) Schizo- 



Lindholm et al.: Schizophrenia-Susceptibility Locus at 6q25 



105 



phrenia: a genome scan targets chromosomes 3p and 8p as 
potential sites of susceptibility genes. Am J Med Genet 60: 
252-260 

Sawcer S, Jones HB, Judge D, Visser F, Compston A, Good- 
fellow PN, Clayton D (1997) Empirical genomewide sig- 
nificance levels established by whole genome simulations. 
Genet Epidemiol 14:223-229 

Sjogren T (1935) Investigations of the heredity of psychoses 
and mental deficiency in two north Swedish parishes. Ann 
Genet 4:253-318 

Schizophrenia Linkage Collaborative Group for Chromosomes 
3, 6 and 8 (1996) Additional support for schizophrenia link- 
age on chromosomes 6 and 8: a multicenter study. Schizo- 
phrenia Linkage Collaborative Group for Chromosomes 3, 
6 and 8. Am J Med Genet 67:580-594 

Straub RE, MacLean CJ, Martin RB, Ma Y, Myakishev MV, 
Harris-Kerr C, Webb BT, O'Neill FA, Walsh D, Kendler KS 
(1998) A schizophrenia locus may be located in region 
10pl5-pll. Am J Med Genet 81:296-301 



Straub RE, MacLean CJ, O'Neill FA, Burke J, Murphy B, Duke 
F, Shinkwin R, Webb BT, Zhang J, Walsh D, Kendler KS 
(1995) A potential vulnerability locus for schizophrenia on 
chromosome 6p24-22: evidence for genetic heterogeneity. 
Nat Genet 11:287-293 

Straub RE, MacLean CJ, O'Neill FA, Walsh D, Kendler KS 
(1997) Support for a possible schizophrenia vulnerability 
locus in region 5q22- 31 in Irish families. Mol Psychiatry 
2:148-155 

Sukumar S, Wang S, Hoang K, Vanchiere CM, England K, 
Fick R, Pagon B, Reddy KS (1999) Subtle overlapping de- 
letions in the terminal region of chromosome 6q24.2-q26: 
three cases studied using FISH. Am J Med Genet 87:17-22 

Turecki G, Rouleau GA, Joober R, Mari J, Morgan K (1997) 
Schizophrenia and chromosome 6p. Am J Med Genet 74: 
195-198 

Zhao JH, Curtis D, Sham PC (2000) Model-free analysis and 
permutation tests for allelic associations. Hum Hered 50: 
133-139 



